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By using a variational approach in combination with the adiabatic approximation we derive a new 
effective ID equation of motion for the axial dynamics of elongated condensates. For condensates with 
vorticity \q\ = or 1, this equation coincides with our previous proposal [A. Munoz Mateo and V. Delgado, 
Phys. Rev. A 77, 013617 (2008)]. We also rederive the nonpolynomial Schrodinger equation (NPSE) 
in terms of the adiabatic approximation. This provides a unified treatment for obtaining the different 
effective equations and allows appreciating clearly the differences and similarities between the various 
proposals. We also obtain an expression for the axial healing length of cigar-shaped condensates and 
show that, in the local density approximation and in units of the axial oscillator length, it coincides with the 
inverse of the condensate axial half-length. From this result it immediately follows the necessary condition 
for the validity of the local density approximation. Finally, we obtain analytical formulas that give the 
frequency of the axial breathing mode with accuracy better than 1%. These formulas can be relevant 
from an experimental point of view since they can be expressed in terms only of the axial half-length and 
remain valid in the crossover between the Thomas-Fermi and the quasi- ID mean-field regimes. We have 
corroborated the validity of our results by numerically solving the full 3D Gross-Pitaevskii equation. 
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I. INTRODUCTION 



In recent years there has been great interest in the physics 
of Bose-Einstein condensates of dilute atomic gases con- 
fined in highly elongated traps fl Q S i & S 0, S H- 
These systems are routinely produced experimentally by 
using microfabricated atom chips [ 10,|TTlll2] or tight opti- 



in terms of an effective equation of lower dimensionality. 
Several pro pos als have been made in recent years in this 
respect IHIH 0, HI SI . 



cal lattices 013111411 and have important applications in gen- 
eration and manipulation of matter-wave solitons 11151 flql 
and in the design of highly sensitive quantum devices such 
as matter-wave interferometers iflvL [lSL U9II . From a the- 
oretical point of view, in the mean-field regime and zero 
temperature limit, they are accurately described in terms 
of a macroscopic wave function 4>(r, t) that satisfies the 
Gross-Pitaevskii equation (GPE) [20] 

iH lti = ("l^ V2 + V{r) + §N m ") ^ (1) 

where A^ is the number of atoms, g = 4:irh 2 a/m is the 
interaction strength determined by the s-wave scattering 
length a, and V(r) is the potential of the confining trap. 
In this work we shall restrict ourselves to condensates with 
repulsive interatomic interactions (a > 0). 

In highly anisotropic cigar-shaped traps the radial con- 
finement can be so tight that the transversal motion be- 
comes practically reduced to zero-point oscillations. Under 
these circumstances, only the slow axial degrees of free- 
dom are relevant and the condensate dynamics becomes ef- 
fectively one-dimensional. This is the quasi- ID mean-field 
regime. More generally, a high anisotropy always induces 
two very different time scales. When the axial motion is 
sufficiently slow in space and time that the radial degrees 
of freedom can adjust practically instantaneously to the dif- 
ferent axial configurations, the radial motion becomes ir- 
relevant and one can still study the condensate dynamics 



By using the adiabatic approximation and the local 
chemical potential that follows from a suitable ansatz for 
the condensate local density [27], we derived in Ref. IJ28I1 
an effective ID equation of motion for the axial dynamics 
of cigar-shaped condensates with repulsive interatomic in- 
teractions. We demonstrated that this equation, which is 
also applicable to condensates containing an axisymmetric 
vortex of topological charge q, is more accurate than pre- 
vious proposals. It also has the advantage that it allows to 
obtain very accurate analytical expressions for a number 
of ground-state properties, and these expressions remain 
valid for condensates with an arbitrary number of particles. 
Despite these merits, from a fundamental point of view it 
would be desirable to find a more systematic way of de- 
riving this equation. In this work, by using a variational 
approach in combination with the adiabatic approximation 
we derive a new effective equation of motion which, for 
\q\ = and 1, coincides exactly with our previous pro- 
posal. We also rederive the nonpolynomial Schrodinger 
equation (NPSE) in terms of the adiabatic approximation. 
This provides a unified method for obtaining the differ- 
ent effective equations and allows appreciating clearly the 
differences and similarities between the various proposals. 
Interestingly, it also demonstrates that in certain cases a 
variational approach based on the chemical-potential func- 
tional can produce more simple and accurate results than 
the usual variational approach based on the energy func- 
tional. We also obtain an expression for the axial heal- 
ing length of cigar-shaped condensates and show that, in 
the local density approximation and in units of the axial 
oscillator length, it coincides with the inverse of the ax- 



ial half-length. From this result it immediately follows the 
necessary condition for the validity of the local density ap- 
proximation. Finally, we obtain analytic formulas that give 
the frequency of the axial breathing mode of an elongated 
condensate with accuracy better than 1% and remain valid 
in the crossover between the Thomas-Fermi (TF) and the 
quasi-lD mean-field regimes. 



II. EFFECTIVE ID EQUATION OF MOTION 



In Ref . 12811 we demonstrated that under usual conditions 
the axial dynamics of highly elongated mean-field conden- 
sates with repulsive interatomic interactions, confined in 
the radial direction by a harmonic potential and contain- 
ing, in general, an axisymmetric vortex of charge q, can be 
described by the effective ID equation 

ih^= ~?t + V z {z)<P + fuv ± J % + iaN |0| V- 



dt 



2m dz 2 



(2) 

This equation incorporates properly the contribution from 
the transverse degrees of freedom through the term propor- 
tional to fvjj_\_. The contribution from the vortex is con- 
tained entirely in the parameter 

3 = - Ug| ' j (3) 

(q (%|)! " {) 

The absence of vortices is a particular case corresponding 
to q = and B q = 1. Since vortices with q > 2 are dy- 
namically unstable 11291 I3QL |3U [3211 . for condensates con- 
taining a multiply quantized vortex the equation above is 
only applicable up to times shorter than the vortex decay 
time. 

Equation (0 was derived in Ref. [28] by applying the 
adiabatic approximation and using for the corresponding 
local chemical potential the analytical expression 



fi_i_(ni) = hw^(\q\ +l-(3 q ) + ^ ± ^//3 g 2 +4ara 1 . (4) 

This expression, in turn, follows from an approximation 
scheme based on a suited TF-like ansatz for the condensate 
local density which, essentially, represents a simple exten- 
sion of the Thomas-Fermi approximation [27]. In what 
follows we will give an alternative derivation that permits 
obtaining the above result in a more systematic way. To 
this end, it is convenient to recall very briefly the main 
steps that led us to Eq. (0. Under usual conditions, the 
time scales characterizing the axial and the radial dynam- 
ics of highly elongated condensates are so different that one 
can appeal to the adiabatic approximation and factorize the 
condensate wave function as l2ll 13311 

ip(r, t) = p(r ± ; m(z, t))(/>(z, t), (5) 

where rj_ = (x,y) and rii(z,t) is the local density per 
unit length along z 

ni (z,t)=N fd 2 r ± ^(r ± ,z,t)\ 2 =N\ ( / ) (z,t)\ 2 . (6) 



Substituting Eq. (0 into the GPE and assuming that the 
axial density varies sufficiently slowly in space and time, 
one obtains [28] 



ih- 



dt 



h 2 d 2 q 
2m dz< 



+ V z (z)<f> + ti X (ni)<f>, (7) 



-^V 2 L + y L (r ± )+ 5 in 1 |v3| 2 j (p = /Uj.(m)(/J, 

(8) 

where we also have assumed a separable confining poten- 
tial V(r) = V ± (r ± ) + V z (z). Equation © shows that the 
axial dynamics is affected by the radial degrees of freedom 
only via the transverse local chemical potential fi±(rii). 
Equation ([8]), which is a stationary GPE, reveals that, at ev- 
ery instant t, the transverse wave function ip coincides lo- 
cally with the equilibrium wave function of an axially uni- 
form condensate characterized by a linear density n\ (z, t). 
This result simply reflects the fact that, for highly elon- 
gated condensates, at every instant of time the (fast) trans- 
verse degrees of freedom can adjust instantaneously to the 
local equilibrium configuration compatible with the axial 
configuration of the condensate. 

In what follows, we shall consider the confining poten- 
tial to be axisymmetric and harmonic in the radial direc- 
tion, while it remains generic in the axial direction 



V(r) 



2 „2 



-mu> i r 



_L'_L 



V z (z). 



(9) 



Using that fJ,x( n i) as given by Eq. © is an accurate ap- 
proximate solution of Eq. ((H), after substituting in Eq. (O 
and taking into account Eq. ©, one finally arrives at Eq. 
©. 

Other different effective equations of motion can be ob- 
tained by using the adiabatic approximation in combination 
with a variational approach. To see this, it is convenient to 
consider the axially uniform condensate described by the 
transverse equation <fS]> as composed of an infinite series of 
identical pieces of length L, with periodic boundary condi- 
tions, and containing N particles each. The exact solutions 
of this equation are the critical points of the energy func- 
tional 

E[<p] 



N 



= d 2 r 



^WM 2 + VxW\ 2 + \gnM 



(10) 
where n\ = N/L is the density per unit length. Thus, the 
problem of finding the eigenfunctions ip that satisfy Eq. |8]l 
is equivalent to the problem of finding, within the whole 
space of admissible functions, those functions that make 
the above energy functional stationary. The corresponding 
local chemical potential then follows from the relationship: 



jUj_(ni) 



dE[p] 
ON ' 



(11) 



In general, however, both problems are equally compli- 
cated, so that, in practice, one usually has to limit the search 



for the critical points of E[<p] to a subspace of convenient 
variational trial functions. The solutions so obtained, in 
general, no longer satisfy exactly the transverse equation 
©, and the corresponding energy E[<p] and chemical po- 
tential fi^(rii) can only be considered as mere estimations 
of the actual values. 

Multiplying by ip* and integrating on the radial coordi- 
nates, Eq. ^ leads to the chemical-potential functional 



/'■± 



\=Jd 2 r ±t p* 



.^-Vi + ^r^ + ffmM'V 

(12) 
Since for condensates with repulsive interatomic interac- 
tions fj,_i [ip] is bounded from below, an independent estima- 
tion for the local chemical potential of the ground state can 
be obtained by minimizing directly the above functional. In 
the ideal-gas perturbative regime (g — ► 0) the system be- 
comes quasi-linear and both estimates coincide. In general, 
however, Eqs. dTUl) and (fT2l lead to different results. When 
the search is performed within the whole space of admis- 
sible functions, minimization of the energy functional al- 
ways yields the correct result. However, as we shall see, 
when the search is restricted to a subspace of variational 
trial functions (as is usually the case), the direct minimiza- 
tion of the functional (PT2l) can lead to a better result for the 
chemical potential of the ground state. 

To the lowest order in the perturbative regime, the 
ground-state solution of Eq. ([8]l compatible with an ax- 
isymmetric vortex of charge q takes the form II2711 



<Pq(r±,0) 



exp(iq9) 

7 



ira ± \q 



(r ± /a ± )""exp(-ri/2ai; 



(13) 

where a± = \/h/muj± is the oscillator length. It is then 
natural to look for the critical points of the energy func- 
tional ( TTOl) within the subspace composed of the above 
functions with the substitution a_L —* ra_i_, where T is a 
dimensionless variational parameter determining the con- 
densate width 04l 13511 - Minimization of Eq. ( TTOb thus 
yields 



1 + 



2an 1 



1/4 



(14) 



Substituting in Eq. ([Tot one obtains the condensate energy 



E , ,, , ,.. L 2an 1 



(15) 



Using this expression in Eq. (fTTT l one finally finds the de- 
sired chemical potential 



3 an i 



m-l = fo*>j-(kl + 1) 



(M + l)/3„ 



1 + 



(16) 



(|g|+i)/3« 



Substitution of Eq. ([TBI into Eq. ([7]) then leads to the fol- 
lowing effective ID equation, to be compared to Eq. (O: 



th- 



at 



h 2 d 2 (f> 

2m dz 2 



+ V z {z)(t>+ 



hu±(\q\ + l) 



1 + 



3aJV|0| 2 
(l?l + l)/3, 



1 + 



2aN\ct>\- 



(17) 



(kl+i)/3, 

This equation, derived here by using the standard adi- 
abatic approximation, is nothing but the nonpolynomial 
Schrodinger equation (NPSE) II36I1 . Equation (TTTl l is also 
applicable to condensates with attractive interatomic inter- 
actions. 

As already said, one can still obtain a different effective 
equation of motion by estimating the chemical potential di- 
rectly from the functional ( PT21) . In doing so, one finds a 
condensate width 

1/4 



1 + 



Acini 



(18) 



and, after substitution in Eq. ( PT21) one arrives at the follow 
ing expression for the local chemical potential: 



u ± = huj ± (\q\ +1),1 + 



Acini 



v vM + i)/V 

Substituting again in Eq. ([7]l one finally obtains 



(19) 



ih m 



h 2 d 2 (j) 

2m dz 2 



+ v z {z)4>+ 



huj ± (\q\ + 1 




4ciN\(j)\' 



(20) 



\q\ + 1)& 

This equation, to be compared to Eqs. (O and (fTTl) . is a 
new effective ID equation governing the axial dynamics 
of the condensate. As can be easily verified, for \q\ = 
and 1 the equation above coincides exactly with our pre- 
vious proposal (0. This is remarkable since the chemical 
potentials that lead to both equations have been derived by 
applying very different techniques: while the chemical po- 
tential ( fT9l ) follows from a variational approach, the chem- 
ical potential (|4]l was derived in Ref. Il27ll by using a suited 
TF-like ansatz for the condensate local density. This non- 
trivial coincidence provides additional support to our pre- 
vious results. Note also that, from a dynamical point of 
view, these cases (corresponding to\q\ =0 and 1) are pre- 
cisely the most relevant ones, since for condensates with 
a larger vortex charge the applicability of the above equa- 
tions is limited to times shorter than the vortex lifetime. 
For \q\ > 2, Eqs. (O and d20l) give somewhat different re- 
sults, which is a consequence of the different way in which 
the chemical potentials (|4| and (fT9l l incorporate the effect 
of the vortex. The point is to determine which of the above 
equations give better results. In Ref. [28] we demonstrated 
that, for q = 0, Eq. Q is somewhat more accurate than Eq. 




FIG. 1: Different theoretical estimates for the local chemical po- 
tential /Zj_ = fi^/hui± as a function of an\. Solid lines have 
been obtained from Eq. @, dashed lines from Eq. 4 19b , and 
dash-dotted lines from Eq. d 1 6b - Open circles are exact results 
obtained by solving numerically Eq.® with no approximations. 



( fTTl l. It remains to be seen whether this is also the case for 
condensates containing a vortex. It is clear that the abil- 
ity of the above equations for reproducing accurately the 
axial dynamics of the condensate is directly related to the 
ability of the corresponding chemical potentials ©, ( fT6l i, 
and ( fT9l l for reproducing accurately the lowest eigenvalue 
of the transverse equation ((8). 

Figured] compares the different theoretical estimates for 
the local chemical potential, obtained from the above for- 
mulas, with exact results (open circles) obtained by solving 
numerically Eq. dU), with no approximations, for a wave 
function of the form exp(iq6)(p(r±) with q = 0, 2 and 4. 
As seen in Ref. [28] our chemical potential, Eq. ©, (solid 
lines) is in good agreement with the numerical results. The 
maximum error is smaller than 1% for q = and 2, and 
smaller than 2.5% for q = 4, and this is so for any value 
of the dimensionless interaction parameter arii (not only 
in the range shown in the figure). As mentioned before, 
the theoretical estimate from Eq. ( fl~9l ) (dashed lines) coin- 
cides exactly with that from Eq. © for q = (and also for 
q = 1, not shown in the figure). For q = 2 the maximum 
error (in the range of the figure) is of the order of 3% and 
for q = 4 it is of the order of 7%. The estimate from Eq. 
( fT6l) (dash-dotted lines) turns out to be somewhat less accu- 
rate. In this case the maximum error is of the order of 5% 
for q = 0, of the order of 7% for q = 2, and of the order of 
10% for q = 4. Moreover, these errors continue increasing 
with arii. These results indicate that the effective ID equa- 
tion (O should give a better description of the condensate 
dynamics than the two other alternative equations. 




FIG. 2: Evolution of {z ) for an elongated condensate with q — 1, 
after a perturbation that excites its axial breathing mode (z = 
z/a z ). Solid lines have been obtained from Eq. (O, dashed lines 
from Eq. fl20| |, and dash-dotted lines from Eq. dl7| >. Open circles 
are exact results obtained from the full 3D GPE. 



To verify this point we have studied the dynamical evo- 
lution of an elongated condensate with A = 0.1 after a 
sudden perturbation that excites its axial breathing mode. 
Figures |2] and [3] show the evolution in time of the mean 
squared axial amplitude (z 2 ) = N^ 1 J dz z 2 rii{z, t) after 
a sudden perturbation at t = in the axial confinement fre- 
quency of the form oo z — > 1.1uj z . Figure|2]corresponds to a 
condensate containing a vortex of charge q = 1 while Fig. 
[5] corresponds to a condensate containing a q = 3 vortex. 
Solid lines have been obtained from our effective equation 
(O, dashed lines have been obtained from the variational 
equation (l20l) . and dash-dotted lines from the NPSE (fT7] l. 
Open circles are exact results obtained by solving numer- 
ically the full 3D GPE with no approximations. Since all 
the three effective equations have the correct perturbative 
(arii <C 1) limit (see inset in Fig. [Ql, they turn out to be 
practically indistinguishable in this regime. As Fig. [2 a) 
shows, for condensates with q = 1 and a (dimensionless) 
peak axial density arii = 1 these equations still give rather 
similar results. The differences increase as arii does, as 
expected from Fig. Q] This can be appreciated from Fig. 
|3b) which corresponds to a condensate with arii = 10. In 
this case Eqs. ([2]) and d20l) . which for q = 1 are indistin- 
guishable, give more accurate results (solid lines) than Eq. 
( fTTl) (dash-dotted lines). 

As Fig. Oa) shows, for q = 3 Eq. (O gives somewhat 




FIG. 3: Evolution of (2 ) for an elongated condensate with q — 3, 
after a perturbation that excites its axial breathing mode (z = 
z/a z ). Solid lines have been obtained from Eq. (O, dashed lines 
from Eq. j20\i , and dash-dotted lines from Eq. dl7| >. Open circles 
are exact results obtained from the full 3D GPE. 



better results than the two other ones even for condensates 
with a peak axial density an\ = 1. In this case, the re- 
sults from Eq. d20l) (not shown in the figure for clarity) lie 
exactly in between the two curves shown. Again, as Fig. 
|3lb) reflects, the differences become more evident as arii 
increases. These results indicate that the effective ID equa- 
tion ((2| is the one that gives a more accurate description of 
the condensate dynamics. It is important to note, however, 
that for condensates with \q\ = or 1 this equation coin- 
cides exactly with the variational equation (l20l l. 

Equation © has the additional advantage that, in com- 
bination with the local density approximation (whose con- 
ditions of validity will be analyzed in detail in Section III), 
it allows to derive a number of useful analytical expres- 
sions for the ground-state properties of elongated conden- 
sates [28]. We next give those that will be needed later on. 
Consider a trapping potential that is also harmonic in the 
axial direction V z {z) = \mu\z 2 , with a z = ^H/mui z 
being the axial oscillator length and A = uj z /ujj_ the trap 
aspect ratio. Then, it can be shown that the total chemical 
potential is given by 12811 



hu>± 



(| g | + l) + i(VA^) 2 , 



(21) 



where Z = Z/a z is the dimensionless axial half-length. 
The axial density profile follows from the formula 



nAz) 



A 



(Vxzy 

4a 



1 



z 2 



(Vxzy 

16a 



1 r 



Z 2 
(22) 
with n\{z) = for \z\ > Z. The interest of the above for- 
mulas lies in the fact that they are valid and accurate (with 
an accuracy typically better than 1%) not only in the TF 
and in the quasi- ID mean-field regime, but also in between 
these two limiting cases. Moreover, since they only depend 
on an easily measurable physical quantity (the condensate 
axial half-length) they can be useful for a more precise ex- 
perimental characterization of these kind of systems. In the 
present work we will use these expressions in order to de- 
rive an approximate analytical formula for the frequencies 
of the axial breathing mode of elongated condensates. 

It can be shown that the axial half-length satisfies the 
following polynomial equation |28|]: 



liVxzr + ^iVxzf 



Xu 



with 



Xi 



XNa 

a± 



(23) 



(24) 



An approximate solution that satisfies the above equation 
for any Xi £ [0, 00), with a residual error [27] smaller than 
0.75% for q = and smaller than 3.2% for 1 < \q\ < 10 
is given by 



y/\Z 



1 



+ 



1 



+ 



1 



(15 Xl )* + ± 57 Xl + 345 (3 X i/&)3 



(25) 
It is clear that the variational equation (l20l also permits 
deriving analytical expressions for the relevant ground- 
state properties. For \q\ = and 1 one obtains the same 
expressions d2T1)-(|25ll. For \q\ > 2, the only modification 
is that the second term on the right-hand side of Eq. (|22] | 
and the first term on the left-hand side of Eq. (l23l become 
multiplied by the factor (3 q /{\q\ + 1). As expected from 
Fig. [Q these modified equations turn out to be somewhat 
less accurate than the above equations, and for this reason, 
except for the present Section, we shall not consider them 
further in this work. Substitution of Eqs. © and (TOl) (with 
a± — > Tax) into Eq. |5]l leads to the equilibrium varia- 
tional wave function 



Y>(r) 



exp(iq0) / r± 
Ta±^/Tr\q\\ V ra -L 



exp(- 



2r 2 ai ■ 



n i( z ) 
N ' 
(26) 



where T(n®) is the (z-dependent) equilibrium condensate 
width given by Eq. ( fT8l) and n\(z) is the axial density pro- 
file given by Eq. d22l) (conveniently modified if \q\ > 2). 
On the other hand, from the equilibrium condensate den- 
sity n(r) = A^|i/;(r)| 2 one finds the following expression 



relating the peak density n(0) with the peak axial density 
n?(0) of a cigar-shaped condensate with no vortices: 



n(0) 



n°(0) 



Tra 2 ±y /l + 4an1(0) 



(27) 



Since an?(0) = (y/X Z/2) 2 + {V\ Z/2) 4 , the above 
formula can also be used to obtain n(0) as a function of 
Z. Equation d27l) has the correct limits in the two ex- 
treme regimes. In particular, in the perturbative regime 
(arii <C 1) it reduces to n(0) = n°(0)/-7ra^, whereas 
in the TF regime (ari\ S> 1) it takes the form n(0) = 
v / n° 1 (0)/a/2Tra 2 L . 



III. AXIAL HEALING LENGTH 

For later purposes, it is convenient to rewrite the effec- 
tive ID equation (O in terms of an equivalent system of 
hydrodynamic equations describing the superfluid dynam- 
ics of the condensate B7I1 . To this end we write the axial 
wave function in polar form 



N(/>(z,t) = y/ ni (z,t)e lS(z ^. 



(28) 



By substituting in Eq. (0, one arrives after some algebra at 
the following equations governing the evolution in time of 
the axial density rii and velocity field v = (h/m)dS/dz: 



drii d . , 



h 2 



(29) 



2 



Im^Jnidz 



dv d ( . . Tr 1 2 

m m + d- z {^ {ni)+Vz+ 2 mv 

(30) 
In these equations, that are completely equivalent to Eq. 
l|2j, the mean-field interaction energy between atoms enters 
through the local chemical potential (ix(rii) given by Eq. 
©. The term proportional to h 2 is the quantum pressure, 
which has its origin in the quantum kinetic energy of the 
system. Note that /J,±(rii) can always be written as 



/ix(ni) = ho±(\q\ + 1) + A_l(™i)- 



(31) 



This can be seen from the transverse equation ©, whose 
exact perturbative solution is an infinite series of the form 
/•*.!_ («i) = ^Ji(kl + 1) + 2f3~ 1 frui ± an 1 + .... From Eq. 
d30l) it is thus clear that the last term on the right-hand side 
of Eq. d3~TT) is the only one that contributes to the dynamics. 
Taking this into account, and defining, as usual, the axial 
healing length £ z as the length scale for which the spatial 
variations in the condensate density induce a quantum pres- 
sure comparable in magnitude to the contribution from the 
interaction energy, i.e., 



h 2 



2mH 



Mx(ni), 



(32) 



one obtains the following expression for the axial healing 
length of elongated condensates: 



e. 



h 



(33) 



y/2m [/zj_(rai) - hw ± (\q\ + 1)] 

In the absence of vortices (q = 0), the above formula re- 
duces to that introduced in Ref. [28]. When the length 
scale A 2 of the spatial variations of the condensate density 
along z is much greater than £ z the contribution from the 
quantum pressure becomes negligible in comparison with 
the contribution from the interaction energy, and Eq. 
reduces to 



dv d ( Tr 1 

rn m + }Tz{ fI±{ni)+Vz+ 2 mv 



0. 



(34) 



In the stationary state v = 0, and the above equation leads 
to 

fi = f i ± (n° 1 ) + V z (z), (35) 

where n\{z) denotes the equilibrium density per unit 
length and ii is the total chemical potential of the conden- 
sate. Equation (l35l l. which holds whenever A z 3> £ z and 
involves no additional approximations, is the so-called lo- 
cal density approximation. 

In what follows, we shall consider a stationary conden- 
sate confined by a trapping potential that is also harmonic 
in the axial direction V z (z) = ^muj 2 z 2 . Substituting Eqs. 
d35l l and (I2TT) into Eq. (I33J I one finds the following ex- 
pression for the local axial healing length of an elongated 
condensate in the local density approximation: 

m = <*/££_. (36 , 

a z y/l - z 2 /Z 2 

This formula shows that at the condensate edges the axial 
healing length diverges, which is an expected result since 
in this region the kinetic energy can never be neglected 
in comparison with the interaction energy. Defining, as is 
usual, the healing length of the condensate as the healing 
length at the center of the density cloud one arrives at 



e. 



z 



(37) 



This formula reflects that in the local density approxima- 
tion and in units of the axial oscillator length, the heal- 
ing length of the condensate coincides with the inverse of 
its axial half-length. As already said, in order for the lo- 
cal density approximation to be valid it is sufficient that 
£ z >C A 2 . Thus, taking into account that for a condensate 
in its stationary state A z is of the order of Z, it follows 
from Eq. ( |37l l that the validity of the local density approx- 
imation requires the condition 

z ( zy 1 

— » — • (38) 

a z \a z J 

It can be easily verified that in the intermediate and TF 
regimes the condition A <C 1 (which always holds for 
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elongated condensates) already guarantees the above re- 
quirement. In the quasi- ID mean-field regime, however, 
the trap aspect ratio A must be sufficiently small so as to 
satisfy the inequality d38l ). Taking into account that in this 
regime an°(0) ~ /3 g (vA Z/2) 2 , one finds that the valid- 
ity of the local density approximation in the perturbative 
(arii ^ 1) regime requires the condition 

4ara?(0) 



A< 



A, 



(39) 



Alternatively, using that for ^ < 1 one has Xi — 
(f3 q /3)(y\Z) 3 , the above condition can also be written 
as 

*«(tr 

IV. COLLECTIVE OSCILLATIONS 

We are interested in obtaining an analytic formula for 
the frequency of the axial breathing mode of elongated 
condensates, applicable also in the crossover between the 
TF and the quasi-lD mean-field regimes. To this end we 
will use the hydrodynamical equation d34b . Considering 
small linear oscillations around the equilibrium configura- 
tion ri\ = n\ + 6rii and v = Sv, after linearization, Eq. 
(f34t takes the form 



d . d ( . n , 



dpi 



dn 1 



5n x + V z 



0, 



Taking into account Eq. (l35t 
to 

d d ( dfi± 

dt dz \ drii 

where, from Eq. ©, one has 

dfi± 2afkv i 



(41) 
the above equation reduces 



Sri! = 0, 



dnx 



^/f3 2 q +4an{ 



(42) 



(43) 



with n\(z) given by Eq. d22l . We look for solutions of 
Eq. (@2J in the form of (axial) dilatations JUIH], that is, 
we introduce a scaling parameter b(t) satisfying the rela- 
tionships z = b(t)z a and v = b(t)z = zb(t)/b(t), and 
assume that the axial density cloud at any t is related to its 
initial value as 

n 1 (,,t) = ^ y n?(, ) = ^ y n?(^ y ). (44) 
For linear oscillations, the case we are interested in, b(t) = 



1 + 5b(t) and 5v = Sb(t)z with 6b(t) ■■ 
the first term in Eq. d42l takes the form 



9 A 

m—ov 

at 



-muj 2 z5b. 



b e~ luJt . Thus, 



(45) 



As for the second term, we note that taking into account 
Eq. (l44l . Srii can be written as 



5rii 



drii 
~~db 



5b 



6=1 



l + Z ^ ! '" {) 



5b. (46) 



Using now Eq. (|22] |. one finds 

(^/XZf 



Srii 



(s/XZ) 2 z 2 
2a Z 2 



\P,+ 




n 



.0 



5b 



From Eqs. @J, (l35l l and d2TT) one has 



z 
'■i ' '■ 2 \ X ~ 'Z 2 

Using this relationship in Eq. d47l ) we obtain 



i 
(47) 

(48) 



5n\ 



{y/X Zf z 2 
2a Z 2 



P 2 + 4an? 



Tli 



5b. (49) 



Substituting now Eqs. ( |43l , d49l and d22l) into the second 
term of Eq. d42l) . after some algebra one finally obtains 

d_ f dfi ± 
dz 



drii 



5n\ ) = m£l 2 (z)uj 2 z5b, 



(50) 



where 



n 2 



2an° 1 (z) 



f3 2 + 4an?(z) 



(51) 



From Eqs. (1451) and ( DOl it is then clear that only in the TF 
and the quasi-lD mean-field regimes Eq. (1421 admit solu- 
tions in the form of dilatations. In the TF regime one has 
4an? 3> p 2 and as a consequence £l 2 (z) = 5/2. In this 
case, Eq. (1421 reflects that the axial breathing mode is an 
oscillatory dilatation with a frequency to 2 = (5/2)uj 2 , in 
good agreement with previous results 113711 . In the quasi- ID 
mean-field regime 4an° <C /3? and one obtains the well- 
known result u; 2 = Z<Jl HUE III]- In between these two 
limiting cases, however, the condensate dynamics does not 
depart too much from a dilatation. This is a direct con- 
sequence of the fact that Q?(z) as given by Eq. d5TT) is a 
slowly varying function of z that can be well approximated 
by its mean value 

-z 






2an { 



dz. 



(52) 



P 2 + 4an° u 

Substituting Eq. (l22l and performing the integration one 
obtains the following analytical expression for the fre- 
quency of the axial breathing mode of a highly elongated 
condensate: 

tanrT 1 (C/V^T^) 



UJ 



ui 



n 2 



+ 



2 2(C 2 + 2) 



+ 



C(C 2 + 2 



3/2 



(53) 

where £ = ^J\j fi q Z. It can be easily verified that the 
above formula has the correct limits in both the TF and the 
quasi-lD mean-field regimes. According to Eq. (l53l one 
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can obtain to 2 experimentally from a measurement of the 
condensate axial half-length Z in the equilibrium configu- 
ration. Alternatively one can also use Eq. ( |25l l to express 
u> 2 /u) 2 z as a function only of Xi an d q - 

One can still derive an independent analytic estimate for 
u 2 by using the formula [6] 

<; 2 } 



" = - 2 W)m- <54> 

This formula, which follows from a sum-rule approach, 
gives an upper bound to the frequency of the axial breath- 
ing mode. Substituting Eq. d22l) in the mean squared am- 
plitude 



(z 2 ) = AT" 1 / dzz 2 n\{z) 



and carrying out the integration, Eq 

u; 2 8(C 2 + 7) 



becomes 



Lo'r 



7-C 2 -(70 + 14C 2 )S' 



(55) 



(56) 



where £ = (X/Z^j dZ/dX. The derivative with respect to 
A can be conveniently rewritten as 



dZ 
~dX 



XVX 



Xi- 



d (V\Z) 
dxi 



l(Vxz) 



(57) 



Now the derivative on the right-hand side can be obtained 
from Eq. (|23T> . In doing so, one finds 

^ (58) 



Xi 

l(3 q (Vxzy 



5Xi 

Substituting this result in Eq. ([56i l and using again Eq. 
to eliminate £ 5 in favor of £ 3 one finally obtains 

4^/ 2 C 3 -15 Xl (C 2 + 5) 



oj 



jj-r 



(59) 

3/3rC 3 -6*i(C 2 + 5) ' 

where £ = \/X/ (3 q Z. As before, taking into account that 
according to Eq. d23l Xi can De rewritten in terms of Z, 
the above equation shows that io 2 can be obtained experi- 
mentally from a measurement of the axial half-length in the 
equilibrium configuration. By using Eq. (I25J I one can also 
rewrite uj 2 /lo 2 as a function only of Xi and f3 q . Equations 
( l53l and Eq. (|59l have been obtained under the assump- 
tion that the quantum pressure has a negligible contribu- 
tion. This implies that in the quasi- ID mean-field regime 
(Xi "C 1) the trap aspect ratio A must be sufficiently small 
so as to satisfy the requirement ( |40l l. In the intermediate 
and TF regimes, however, the condition A < 1 already 
guarantees the inequality ( |38l l to hold true. It can be eas- 
ily verified that Eqs. (l53l l and ( |59l coincide each other 
within 0.65% for q = and within 0.77% for q = 1. 
In order to determine to what extent the above formulas 
reproduce the experimental results we have performed a 
computer experiment based on the numerical solution of 
the full 3D Gross-Pitaevskii equation ([]]). We have con- 
sidered 87 Rb condensates at zero temperature confined in a 
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FIG. 4: Evolution of a Rb condensate with \i — 10 and no 
vortices, in a trap with A = 0.1 after a 10% fluctuation in the 
axial trap frequency. Top: isosurfaces of the condensate density 
taken at intervals At = Q.Slu^ 1 . Bottom: evolution of the mean 
squared axial amplitude (z = z/a z ). Open circles correspond to 
the numerical data and the solid line is the best sinusoidal fit to 
these data. 



magnetic trap with axial frequency oo z = 2n x 3 Hz and 
different radial frequencies uj± = co z /A with A < 1/10. 
These condensates have a scattering length a = 5.29 nm 
and an oscillator length a z = 6.23 /im. For condensates 
with a given vortex charge q, we vary the relevant parame- 
ter Xi = yXNa/a z by modifying the trap aspect ratio A 
and/or the number of atoms N. For a given value of Xi we 
always start with A = 1/10 and calculate the condensate 
equilibrium configuration in the trap. We then excite the 
axial breathing mode by introducing a small fluctuation in 
the axial frequency for a period of time t = lo^ 1 and let 
the condensate oscillate in the trap. At this stage we moni- 
tor the mean squared axial amplitude as a function of time 
and extract the corresponding oscillation frequency from 
a nonlinear least-squares fit of a sinusoidal function to the 
numerical data. We have considered different perturbations 
of the trap frequency, ranging from 1% to 10%, and differ- 
ent fitting intervals and have always obtained the same re- 
sults within our required precision (better than 0.1%). Fig- 
ure |4] displays an example corresponding to a condensate 
with 37220 atoms (xi = 10) and no vortices, in a trap 
with an aspect ratio A = 1/10. The top panel shows the 
condensate evolution after a 10% perturbation in the axial 
trap frequency. The images are isosurfaces of the conden- 
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FIG. 5: Squared frequency of the axial breathing mode (w/wz) 
of highly elongated condensates with q = and 1, as a function 
of xi- Open symbols are numerical results obtained from the full 
3D GPE. Solid curves are the theoretical prediction from Eq. 
and dashed curves are the theoretical prediction from Eq. 



sate density (corresponding to 5% of the maximum den- 
sity) taken at intervals At = O.Su;^ 1 . The bottom panel 
shows the evolution of the mean squared axial amplitude in 
units of a 2 z (z stands for z/a z ). Open circles correspond to 
the numerical data and the solid line is the best sinusoidal 
fit to these data. From this fit we obtain an oscillation fre- 
quency ijj = 1.605u; 2 with an accuracy better than 0.1%. 

In general, a non-negligible contribution of the quantum 
pressure manifests itself as an appreciable dependence of 
the oscillation frequency on the parameter A. For a given 
value of Xi, as A increases the contribution of the quantum 
pressure also increases, producing in turn an increase in 
the condensate oscillation frequency. This fact can be used 
in particular to quantify the contribution of the quantum 
pressure or, equivalently, the validity of the local density 
approximation. To guarantee that the numerical results ob- 
tained are comparable with the above analytic formulas we 
have repeated the procedure explained before using smaller 
and smaller values of A until reaching an oscillation fre- 
quency exhibiting no appreciable dependence on A. An 
estimation of the values required can be obtained from Eq. 
d40l) . Figure [5] shows the (squared) frequency of the axial 
breathing mode (lu/uj z ) 2 of highly elongated condensates 
with vorticity q = and 1, as a function of Xi- Open 
symbols are numerical results obtained from the full 3D 
GPE with an accuracy of the order of 0.1%. Solid curves 
are the theoretical prediction from Eqs. (l53l and d25l) and 
dashed curves are the theoretical prediction from Eqs. ( |59l l 
and (l25l l. Equation (l53l reproduces the numerical results 
with an accuracy better than 0.8% for q = and better 
than 0.9% for q = 1. Equation ( |59l l is somewhat more 
accurate. It reproduces the numerical results with an ac- 
curacy better than 0.4% for q = and better than 0.65% 



for q = 1. As Fig. [5]reflects, in the crossover between the 
TF and the quasi- ID mean-field regimes there exists a clear 
dependence of the frequency to on the vortex charge q. This 
fact could be used for an indirect detection of unit-charge 
vortices in the crossover regime. Because of their accuracy 
and range of applicability the above formulas could also be 
useful from an experimental point of view, for instance, for 
calibration of the trap frequencies. 



V. CONCLUSION 

In this work, by using the adiabatic approximation in 
combination with a variational approach for determining 
the local chemical potential associated with the corre- 
sponding radial equation, we have derived different ef- 
fective ID equations of motion for the axial dynamics of 
highly elongated condensates. We have shown that the 
minimization of the radial energy functional within a sub- 
space of convenient variational trial functions leads to the 
NPSE. We also have demonstrated that in certain cases a 
variational approach based on the chemical-potential func- 
tional can produce better results than the usual variational 
approach based on the energy functional. In this regard, 
while the minimization of the energy functional within 
the whole space of admissible functions always yields the 
correct chemical potential, the same does not remain true 
when this minimization is restricted to a certain variational 
subspace. In fact, we have seen that in this case it is possi- 
ble to obtain a more simple and yet more accurate result by 
minimizing directly the chemical-potential functional. By 
doing this, we have derived a new effective ID equation of 
motion which, for condensates with vorticity \q\ = and 
1, coincides exactly with our previous proposal, obtained in 
Ref. [28] by using a suited TF-like ansatz. The variational 
approach followed in this work provides us with a unified 
method for obtaining the different effective ID equations 
in a systematic way and permits us to appreciate clearly the 
differences and similarities between the various proposals. 
A direct comparison with numerical results from the full 
3D GPE indicates that the effective ID equation proposed 
in Ref. [28] is the one that gives a more accurate descrip- 
tion of the condensate axial dynamics. This equation also 
has the advantage that, in combination with the local den- 
sity approximation, it allows to derive accurate analytical 
formulas for a number of ground-state properties. These 
formulas are valid and accurate (with an accuracy typically 
better than 1%) not only in the TF and in the quasi- ID 
mean-field regime, but also in between these two limiting 
cases. Moreover, since they only depend on an easily mea- 
surable physical quantity (the condensate axial half-length) 
they can be useful for a more precise experimental charac- 
terization of these kind of systems. 

We also have obtained an expression for the axial heal- 
ing length of elongated condensates and have found that, 
in the local density approximation and in units of the axial 
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oscillator length, it coincides with the inverse of the con- 
densate axial half-length. From this result it immediately 
follows the necessary condition for the validity of the local 
density approximation. 

Finally, we have obtained approximate analytic formu- 
las that give the frequency of the axial breathing mode of 
an elongated condensate with accuracy better than 1% and 
remain valid in the crossover between the TF and the quasi- 
1D mean-field regimes. These formulas, which as the rest 
of the ground-state properties can be expressed in terms 
only of the axial half-length Z, could be relevant from an 
experimental point of view since in the crossover regime 
(the regime where most experiments are carried out) the 
usual formulas are not applicable. 
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